Global Phase Diagram for Magnetism and Lattice Distortion of 

Fe-pnictide Materials 
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We study the global phase diagram of magnetic orders and lattice structure in the Fe-pnictide 
materials at zero temperature within one unified theory, tuned by both doping and pressure. On 
the low doping and high pressure side of the phase diagram, there is one single transition, which is 
described by a z = 2 mean field theory with very weak run-away flows; on the high doping and low 
pressure side the transition is expected to split to two transitions, with one 0(3) spin density wave 
transition followed by a z = 3 quantum Ising transition at larger doping. The fluctuation of the 
strain field fluctuation of the lattice will not affect the spin density wave transition, but will likely 
drive the Ising nematic order transition more mean field like through a linear coupling, as observed 
experimentally in BaFe2- x Co x As2. 



I, INTRODUCTION 

The Iron-superconductor, for its potential to shed new 
light on the non-BCS type of superconductors, has at- 
tracted enormous interests since early this year. Despite 
the complexities and controversies on the superconduct- 
ing mechanism, the minimal tight-binding model, or even 
the exact pairing symmetry of the cooper pair, these 
samples do share two common facts: the tetragonal- 
orthorhombic lattice distortion and the (n, 0) spin den- 
sity wave (SDW) [l|. Both effects are suppressed under 
doping and pressure, and they seem always to track each 
other in the phase diagram. In Ref . [ljl [l6| , the lattice 
distortion is attributed to preformed spatially anisotropic 
spin correlation between electrons, without developing 
long range SDW i.e. the lattice distortion and SDW both 
stem from magnetic interactions. More specifically, the 
Ising order parameter a is represented as a — <j>\ ■ <p2 , 4>i 
and </>2 are two Neel orders on the two different sublat- 
tices of the square lattice. 

Since this order deforms the electron Fermi surface, 
equivalently, it can also be interpreted as electronic ne- 
matic order. The intimate relation between the structure 
distortion and SDW phase has gained many supports 
from recent experiments. It is suggested by detailed X- 
ray, neutron and Mossbauer spectroscopy studies that 
both the lattice distortion transition and the SDW tran- 
sition of LaFeAs(Oi_ x F x ) are second order [l7|, where 
the two transitions occur separately. However, in un- 
doped AFe2As2 with A = Sr, Eu, Ba, Ca, the structure 
distortion and SDW occur at the same temperature, and 
the transition becomes a strong first order transition 
[H, 0, H, 0, 0]- Also, recent neutron scattering measure- 
ments on Fei+ySe x Tei_ x indicate that in this material 
the SDW wave vector is (tt/2,tt/2) for both sublattices 
24j instead of (tt, 0) as in 1111 and 122 materials, and the 
low temperature lattice structure is monoclinic instead of 
orthorhombic (choosing one-Fe unit cell). These results 
suggest that the SDW and structure distortion are indeed 
strongly interacting with each other, and probably have 



the same origin. The sensitivity of the location of the 
lattice distortion transition close to the quantum critical 
point against the external magnetic field (magnetoelastic 
effect) can further confirm this unified picture. 

The clear difference between the phase diagrams of 
1111 and 122 materials can be naturally understood in 



We can write 



the unified theory proposed in Ref. 
down a general Ginzburg-Landau mean field theory for 
a, 4>\ and (p2 ■ 



Fgl = (V m ct) 2 +r a a 2 + ^( V ^) 2 +r^\$ a \ 



uacpx ■ <p2 + ■ ■ ■ 



(1) 



r a and are tuned by the temperature. For a purely 
two dimensional system, the Ising order which induces 
the lattice distortion occurs at a temperaure controlled 
by the inplane spin coupling Ti S i ng ~ Jj„, while there is 
no SDW transition at finite temperature; for weakly cou- 
pled two dimensional layers, the Ising transition temper- 
ature is still of the order of the inplane coupling, while the 
SDW transition temperature is T s( i w ~ ■/in/ln(Jj n /J z ), 
with J z <C Ji n representing the intcrlayer coupling. This 
implies that on quasi two dimensional lattices, Ti S i ng > 
T s dw, or Ar = — r a is large in the Ginzburg-Landau 
mean field free energy in Eq. [TJ In real systems, The 1111 
materials are much more anisotropic compared with the 
122 materials, since the electron band structure calcu- 
lated from LDA shows a much weaker z direction disper- 
sion compared with the 122 samples [26j]; also the upper 
critical field H C 2 of 122 samples is much more isotropic 
[25l | . This justifies treating the 1111 materials as a quasi 
two dimensional system, while treating the 122 materi- 
als as a three dimensional one. When J z and Ji„ are 
close enough, Ar is small, and the interaction between 
the Ising order parameter and the SDW will drive the 
transition first order by minimizing the free energy Eq. 
[U The phase diagram of free energy Eq. Q] is shown in 
Fig.ffl 

Motivated by more and more evidences of quan- 




FIG. 1: The schematic phase diagram of Ginzburg-Landau 
mean field theory in Eq. [TJ plotted against r — r a + r#, and 
Ar = — r a . r is linear with temperature T, while Ar is 
tuned by anisotropy ratio J z /Ji n - When Ar is small, the 
interaction between cj>\ and <j>2 induces a strong first order 
transition, which corresponds to the undoped 122 materials 
with more isotropic electron kinetics; when Ar is large, the 
transition is split into two transitions, with an Ising transition 
followed by an SDW transition at lower temperature, and this 
is the case in the 1111 materials with quasi two dimensional 
dispersions. The multicritical point Ar c is determined by u. 



turn critical points in the Fe-pnictides superconductors 

@, B H, EH IH El, in this work > we wil1 explore the 
global phase diagram of magnetic and nematic orders at 
zero temperature, tuned by two parameters, pressure and 
doping. In section II we will study the phase diagram for 
quasi two dimensional lattices, with applications for 1111 
materials, and in section III the gear will be switched to 
the more isotropic 3d lattices of 122 materials. Section IV 
will briefly discuss the effect of the coupling between the 
Ising transition and the strain tensor of the lattice, which 
will drive the finite temperature Ising nematic transition 
a mean field transition, while the SDW transition re- 
mains unaffected, as observed in BaFe2_ x Co x As2 [Tol ] . 
The analysis in our current work are all only based on 
the symmetry of the system, and hence independent of 
the details of the microscopic model. 



II, QUASI TWO DIMENSIONAL LATTICE 

In Ref . [15J , the zero temperature quantum phase tran- 
sition was studied for weakly coupled 2d layers with finite 
doping. Since the hole pockets and the electron pockets 
have small and almost equal size, the slight electron dop- 
ing would change the relative size of the electron and 
hole pockets substantially. Also, the neutron scattering 
measurement suggests that the SDW order wave vector 
is independent of doping in 1111 materials [2|. Therefore 
under doping the low energy particle-hole pair excitations 
at wave vector (tt, 0) are lost very rapidly, and the spin 
density wave order parameter at low frequency and long 
wavelength limit can no longer decay with particle-hole 















FIG. 2: a, The schematic two dimensional fermi pockets of 
1111 materials, the red circles are two very close hole pockets, 
the blue ovals are electron pockets, b — d, the relative posi- 
tion of hole and electron pockets after translating by (tv, 0) 
and (0, 7r) in the momentum space, in the low doping, critical 
doping, and high doping regimes. 



pairs (the fermi pockets are schematically showed in Fig. 
[2]). After integrating out electrons we would obtain the 
following z — 1 Lagrangian (l5j : 

2 

L = dJt-dJt-rfi+u^ + L', 

t=l n=T,x,y 

V = 1 $ 1 (€-9l)-$2+lSi\ 2 \h\ 2 -a($ 1 -$ 2 ) 2 {2) 



which contains no damping term. The first three terms 
of the Lagrangian describe the two copies of 3D 0(3) 
Neel orders on the two sublattices. The a term is the 
only relevant term at the 3D 0(3) transition, since it has 
positive scaling dimension A [a] = 0.581 We expect 
this term to split the two coinciding 0(3) transitions into 
two transitions, an 0(3) transition and an Ising transition 
for Ising variable a = <pi ■ (f>2 , as observed experimentally 
in 1111 materials Q. 

The two transitions after splitting are an 0(3) transi- 
tion and an Ising transition. The 0(3) transition belongs 
to the 3D 0(3) universality class, while the Ising tran- 
sition isaz = 3,e? = 2 mean field transition. This is 
because the Ising order parameter does not double the 
unit cell, and hence can decay into particle-hole pairs at 
momentum (0,0). The standard Hertz-Millis theory [34[ 
would lead to a z = 3 mean field transition [3, [3l[ . 

Now let us turn on another axis in the phase diagram: 
the pressure. Under pressure, the relative size of hole and 
electron pockets are not expected to change. Therefore 
under translation of Q = (w, 0) in the momentum space, 
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FIG. 3: The plot of u in Eq. [8] against anisotropic dispersion 
coefficient 7, between the isotropic limit 7 = to anisotropic 
dispersion with 7 = 1.95 



the hole pocket will intersect with the electron pocket 
(Fig. which leads to overdamping of the order pa- 
rameters. The decay rate can be calculated using Fermi's 
Golden rule: 



Im[x(w,<?)] 



d 2 k 



[f(e k+q ) - f(e k+<3 )] 



(2tt) 2 w J ^' k + 

x 5{u) - e k+q + e fc+ g)|(fc + Q\$i,q\k + q)\ 



CO, 



\v h x v e 



(3) 



Vh and t> e are the fermi velocity at the points on the hole 
and electron pockets which are connected by wave vector 
(it, 0). The standard Hertz-Millis 34j formalism leads 
to a coupled z — 2 theory in the Euclidean momentum 
space with Lagrangian 

2 

L i = ^^•(|w| + 9 2 + r)0' l +70 1 ( (? 2_ (Z 2 ) .0 2 + L / ! 

i=l 

L' = A(\$ 1 \ 4 + \$ 2 \ 4 )-a($ 1 -$ 2 ) 2 + C\$ 1 \ 2 \$ 2 \ 2 . (4) 

The parameter r can be tuned by the pressure. The Ising 
symmetry of a = <pi ■ 4> 2 on this system corresponds to 
transformation 



y, y^x, 



a — > —a. 



(5) 



This Ising symmetry forbids the existence of term <f>\ ■ <f> 2 
in the Lagrangian, while the mixing term j4>i(q 2 — q 2 ) -4> 2 
is allowed. 

We can diagonalize the quadratic part of this La- 
grangian by defining <\>a — {4>i + <p 2 )/^/2 and <pB = 
{$! - $ 2 )/V2: 



L q = 



^•(M + (l-^ + (l + ^ + r)^4 



2 )q y 



L', 



L' = A(\$ A \ 4 + \$ B \ 4 ) + B($ A -$ B ) 2 + C\$ A \ 2 \$ B \$) 
After the redefinition, the Ising transformation becomes 



2/, y^x, 



(7) 



Naively all three quartic terms A, B and C are marginal 
perturbations on the z — 2 mean field theory, a coupled 
renormalization group (RG) equation is required to de- 
termine the ultimate fate of these terms. Notice that 
the anisotropy of the dispersion of (f>A and <f>B cannot 
be eliminated by redefining space and time, therefore the 
number 7 will enter the RG equation as a coefficient. The 
final coupled RG equation at the quadratic order for A, 
B and C reads: 



dA 
dln7 
dB 

din7 

dC 
dhd 



-22A 2 - -B 2 - -C 2 - BC. 
2 2 



-5uB 2 - SAB - 8uBC, 



-uB 2 - AAB - 20 AC - 4uC 2 . 



(8) 



u is a smooth function of 7, which decreases smoothly 
from u = 1 in the isotropic limit with 7 = to u = in 
the anisotropic limit with 7 = 2 (Fig J3]). The self-energy 
correction from the quartic terms will lead to the flow 
of the anisotropy ratio 7 under RG, but the correction 
of this flow to the RG equation Eq. [5] is at even higher 
order. 

The typical solution of the RG equation Eq. [H] is plot- 
ted in Fig. |U One can see that the three parameters A, 
B and C all have run-away flows and eventually become 
nonperturbative, and likely drive the transition weakly 
first order. However, the three coefficients will first de- 
crease and then increase under RG flow. This behav- 
ior implies that this run-away flow is extremely weak, 
or more precisely even weaker than marginally relevant 
perturbations, because marginally relevant operators will 
still monotonically increase under RG flow, although in- 
creases slowly. Therefore in order to see this run-away 
flow, the correlation length has to be extremely long i.e. 
the system has to be very close to the transition, so 
the transition remains one single second order mean field 
transition for a very large length and energy range. At 
the finite temperature quantum critical regime, the stan- 
dard scaling arguments lead to the following scaling laws 
of physical quantities like specific heat and the spin lat- 
tice relaxation rate of NMR contributed by the quantum 
critical modes 14|: 



a~Tln(i), i~ Const. 



(9) 



These scaling behaviors are obtained from ignoring the 
quartic perturbations. The quartic terms are marginal 
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FIG. 4: The solution of the RG equation Eq. IS All three 
quartic perturbations decrease first, then increase and finally 
become nonperturbative. The run-away flow is weaker than 
marginally relevant perturbations. 



for a rather large energy scale (Fig. therefore to pre- 
cisely calculate the physical quantities one should per- 
form a perturbation theory with constant A, B and C, 
which may lead to further logarithmic corrections to the 
scaling laws. 

The z direction tunnelling of 4>a and 4>b between layers 
has so far been ignored, which is also a relevant perturba- 
tion at the z = 2 mean field fixed point. The z direction 
tunnelling is written as J z (f> a . z • 4>a, z +i, which has scaling 
dimension 2 at the z = 2, d = 2 mean field fixed point, 
and it becomes nonperturbative when 



Jin ^ (— ) 2 
J z a 



(10) 



This equation implies that if the tuning parameter r is 
in the small window r < J z /Ji n , the transition crossover 
back to a z = 2, d = 3 transition, where all the quartic 
perturbation A, B and C are irrelevant. Since at the two 
dimensional theory these quartic terms are only weakly 
relevant up to very long length scale, in the end the in- 
terlayer coupling J z may win the race of the RG flow, 
and this transition becomes one stable mean field second 
order transition. 

Now we have a global two dimensional phase diagram 
whose two axes are doping and pressure. The two second 
order transition lines in the large doping and low pressure 
side will merge to one single mean field transition line 
in the low doping and high pressure side of the phase 
diagram. Then inevitably there is a multicritical point 
where three lines merge together. At this multicritical 
point, the hole pockets will just touch the electron pocket 
after translating by wave vector {%, 0) (Fig. ^Bp)- Now 
the SDW order parameter cf>A and 4>b can still decay 
into particle-hole pairs, the Fermi's Golden rule and the 
lattice symmetry lead to the following overdamping term 
in the Lagrangian: 



L q = ( 




Disorder 



Nematic 




FIG. 5: The global phase diagram of quasi two dimensional 
materials, with applications for 1111 materials. The finite 
temperature transition is always split to an Ising nematic 
transition and a SDW transition. The zero temperature tran- 
sitions depend on the doping and pressure. In the high dop- 
ing and low pressure side, the transition is split to two, as 
observed in experiments; in the low doping and high pressure 
side, there is one single transition very close to the mean field 
solution. A multicritical point where the three transition lines 
merge is identified, which is expected to be a strongly coupled 
fixed point. 



g is a constant, which is in general not unity because 
the system only enjoys the symmetry Eq. [7] The naive 
power-counting shows that this field theory has dynami- 
cal exponent z = 5/2, which makes all the quartic terms 
irrelevant. However, since the hole pockets and electron 
pockets are tangential after translating (n, 0), the ex- 
pansion of the mean field free energy in terms of the 
order parameter cf>A and 4>b contains a singular term 
L s ~ 10,4 1 5 / 2 + |</>_b| 5 / 2 , which becomes very relevant at 
this naive z = 5/2 fixed point. Similar singular term was 
found in the context of electronic nematic-smectic tran- 
The existence of this singular term implies 



21|. 



)\$b? 



(11) 



sition 

that, it is inadequate to start with a pure Bose theory 
by integrating out fermions, one should start with the 
Bose-Fermi mixed theory, with which perform the RG 
calculation. We will leave this sophisticated RG calcula- 
tion to the future work, right now we assume this multi- 
critical point is a special strongly interacting fixed point. 
The schematic three dimensional global phase diagram is 
shown in Fig. [5] 

In real system, due to the more complicated shape 
of the electron and hole pockets, with increasing doping 
the pockets will experience cutting and touching several 
times after translating by (it, 0) in the momentum space. 
We have used a five-band model developed in Ref. [1(3] 
with all the d— orbitals on the Fe atoms, and calculated 
the mean field phase diagram close to the critical dop- 
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FIG. 6: Numerical results of r c of 4> a due to coupling to elec- 
trons, x axis is the electron doping. The peak of this curve 
corresponds to the critical doping x c — 7.6% where electron 
and hole pockets touch each other after translating the hole 
pockets by the SDW wave vector. The two pockets intersect 
(separate) if doping is smaller (larger) than this critical dop- 
ing. If a; > x c , the transition is split into two transitions by 
quantum fluctuations; if x < x c , the transition is a z = 2, 
d — 2 transition with a very weak run-away flow in 2d. 
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FIG. 7: The plot of the hole and electron pockets after trans- 
lating the hole pockets by the SDW wave vector, at the crit- 
ical doping x — 7.6%. The green circle is electron pocket 
located around (0, n) and the blue one is electron pocket lo- 
cated around (7r,0). 



ing. The order parameter <f> a couple to the electrons at 
the Fermi surface as: J2k $ a ' C k^ c k+Q H- c - The mean 
field energy of electrons due to nonzero spin order pa- 
rameter <j> a will renormalizc r in field theory Eq. [6l and 
hence the critical r c depends on the shape of the Fermi 
surface, which is tuned by doping. The critical r c is ex- 
pected to be proportional to the critical pressure p c in 
the global phase diagram. r c as a function of doping is 
plotted in Fig. [5J and the shapes of the Fermi pockets at 
the critical doping x — 7.6% is plotted in Fig. [7J 

The z = 2 quantum critical behavior discussed in this 
section is only applicable to small enough energy scale. 
First of all, the damping term \uj\ always competes with 
a quadratic term to 2 in the Lagrangian, and at small 
enough energy scale the linear term dominates. If we 
assume the coupling between the spin order parameter 
(f> a and the electrons is of the same order as the effec- 
tive spin interaction J, the damping rate is linear with 
~ J 2 ui/E 2 , while the quadratic term is ~ lo 2 / J . There- 
fore the frequency should be smaller than J 3 /E 2 in order 
to apply the z — 2 field theory Eq. [S) The value of J 
has been calculated by LDA 11811 . and also measured with 
inelastic neutron scattering [19j , and both approaches in- 
dicate that J ~ 50meV. Ef is the Fermi energy of the 
Fermi pockets, which is of the order of 200meV. There- 
fore the frequency-linear damping term will dominate the 
the frequency quadratic term in the Lagrangian as long 
as ui < 3meV. 

The damping rate of order parameters 4> a is calculated 
assuming the Fermi surface can be linearly expanded 
close to the intersection point after translation in the 
momentum space, the criterion to apply this assumption 



depends on the details on the Fermi surface. In the par- 
ticular situation under discussion, this crossover energy 
scale is to ~ 30meV in the undoped material, which is 
larger than the upper limit of 3meV we obtained previ- 
ously. Therefore the ultraviolet cut-off of field theory Eq. 
[5] is estimated to be 3meV. 

Ill, THREE DIMENSIONAL LATTICES 

As mentioned in the introduction, compared with 
the 1111 materials, the 122 materials are much more 
isotropic, so we will treat this family of materials as 
a three dimensional problem. If after translation by 
(tt, 0) the hole pockets intersect with the electron pock- 
ets, the zero temperature quantum transition is described 
a z = 2, d = 3 transition with analogous Lagrangian as 
Eq. [U which becomes a stable mean field transition. The 
finite temperature transition is described by two copies 
of coupled 3D 0(3) transition. If the finite tempera- 
ture transition is split into two transitions close to the 
quantum critical point, as observed in BaFe2- x Co x As2 
|ld |. one can estimate the size of the splitting close to 
the quantum critical point. These two transitions, as 
explained before, is driven by the only relevant pertur- 
bation a(<f>i ■ 02) 2 at the coupled 3D 0(3) transition, be- 
cause the Ising order parameter is obtained by minimiz- 
ing this term through Hubbard-Stratonovich transforma- 
tion. The scaling dimension of a at the 3D 0(3) transi- 
tion is A [a] = 0.581, while a at the z = 2, d = 3 mean 
field fixed point has dimension —1. Therefore close to the 
quantum critical regime, to estimate the effect of a one 



6 



should use the renormalized value ccr ~ o^ -1 ~ ar 1 / 2 . 
The size of the splitting of the finite temperature transi- 
tion close to the quantum critical point can be estimated 
as 



1/("AH) 



Q ,l/( I /A[c.]) r l/(2 I /A[c.]) < f 12 \ 



v is the exponent defined as £ ~ t~ u at the 3D 0(3) 
universality class. T c still scales with r in terms of a 
universal law T c ~ r z/(<i-2+z)^ fp ne number a can be 
estimated in a Ji — J2 Heisenberg model on the square 
lattice as introduced in Ref. (HI, the value is given by 
a ~ J\jJ\- However, Ji — J 2 model is not designed for 
describing a metallic phase, so the legitimacy of applying 
the J\ — J2 model to Fe-pnictides is still under debate. 
In the finite temperature quantum critical regime, the 
specific heat, NMR relaxation rate 1/T\ scale as 



J-l 



(13) 



The similar analysis also applies when the finite tem- 
perature transition is one single first order transition, 
which is the more common situation in 122 materials. 
One can estimate the jump of the lattice constant, and 
the jump of the SDW order parameter at the finite tem- 
perature first order transition close to the quantum crit- 
ical point as 



8 a 



a /3/(uA[a]) r 0/{2uA[a]) 



a 2f}/{ v A[a]) r P/(vA[ a ))^ 



(14) 



a is the lattice constant, which is linear with the Ising 
order parameter 4>\ ■ 4>2- /? is the critical exponent at the 
3D 0(3) transition defined as {4> s dw) ~ ^ ■ 

If under doping the hole pockets and electron pockets 
do not intersect (which depends on the details of z direc- 
tion dispersions), this transition becomes two copies of 
coupled z = 1, d = 3 transition with three quartic terms 
A, B and C: 



L„ 



7 U 2 



">S„2 



(w ' + + M + (1-77K + qt+r 



L' 



L' 



A(\$a\ 4 



B($ A ■ $ B ) 2 + C\$ A \ 2 \Uf5) 



The coupled RG equation of the quartic terms are exactly 
the same as the one in Eq. [5J therefore this free energy 
is also subjected to an extremely weak run-away flow, 
which is negligible unless the length scale is large enough. 
Again, one can estimate the universal scaling behavior in 
the quantum critical regime contributed by the quantum 
critical modes: 



1 

T. 



(16) 



IV, COUPLING TO A SOFT LATTICE 

Recent specific heat measurement on BaFe2- x Co x As2 
reveals two close but separate transitions at finite tem- 
perature, with a sharp peak at the SDW transition, and 
a discontinuity at the lattice distortion transition [loj ]. 
A discontinuity of specific heat is a signature of mean 
field transition, in contrast to the sharp peak of Wilson- 
Fisher fixed point in 3 dimensional space. The specific 
heat data suggest that the nature of the Ising nematic 
transition is strongly modified from the Wilson-Fisher 
fixed point, while SDW transition is unaffected. In the 
following we will attribute this difference to the lattice 
strain field fluctuations. 

The SDW transition at finite temperature should be- 
long to the 3D 0(3) transition ignoring the lattice. The 
0(3) order parameter <j> couples to the the lattice strain 
field with a quadratic term [33| : 



\<t>\ 2 (d x u x + dyUy + X'd z u z ), 



(17) 



which after integrating out the displacement vector gen- 
erates a singular long range interaction between \<f>\ 2 in 
the real space: 



(18) 



/ is a dimensionless function which depends on the direc- 
tion of f—f". The scaling dimension of g is A[g] = 2/v—3, 
and v is the standard exponent at the 3D 0(3) transition, 
which is greater than 2/3 according to various types of 
numerical computations [32I ]. Therefore this long range 
interaction is irrelevant at the 3D 0(3) transition, and by 
coupling to the strain field of the lattice, the SDW tran- 
sition is unaffected. However, if the SDW has an Ising 
uniaxial anisotropy, the SDW transition becomes a 3D 
Ising transition with v < 2/3, and the strain field would 
lead to a relevant long range interaction. 

However, since the symmetry of the Ising order param- 
eter a is the same as the shear strain of the lattice, the 
strain tensor will couple to the coarse grained Ising field 
<I> as 



= \$(d x Uy + d y u x ) 



(19) 



u is the displacement vector, The ellipses are all the elas- 
tic modulus terms. Notice that we have rotated the co- 
ordinates by 45 degree, since the true unit cell of the 
system is a two Iron unit cell. After integrating out the 
displacement vector u, the effective free energy of $ gains 
a new singular term at small momentum: 



(20) 



/ is a function of spherical coordinates 9 and (j> defined as 
(k x ,k y , k z ) = fc(cos(#) cos(0), cos(0) sin(^), sin(0)), but / 
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is independent of the magnitude of momentum k. By 
tuning the uniform susceptibility r, at some spherical 
angle of the space the minima of / start to condense, 
we will call these minima as nodal points. These nodal 
points are isolated from each other on the two dimen- 
sion unit sphere labelled by the solid angles 8, <fi, and 
are distributed symmetrically on the unit sphere (8, cj>) 
according to the lattice symmetry transformation. Now 
suppose one nodal point of / is located at (8 , (j) ),we ro- 
tate the z direction along (#o,</>o), and expand / at this 
nodal point in terms of 8 = 8 — 8q, the whole free energy 
can be written as 



as $2 and $2, the expanded free energy reads 



F= / q 2 dq8d8(q 2 + X8 2 + r)|* 



0($ 4 ). (21) 



Notice that if f(8o, <fro) is a nodal point, then /(tt—Oq, ir+ 
<j>o) has to be another nodal point. The naive power 
counting shows that effectively the spatial dimension of 
this field theory Eq. [5l]is D — 5, and the scaling dimen- 
sion of Q qt g is —7/2. The quartic term <i> 4 takes an un- 
usual form in the new momentum space of q, 8, but the 
straightforward power counting indicates that it is still 
an irrelevant operator. Therefore the strain tensor fluc- 
tuation effectively increases the dimension by two, which 
drives the transition a mean field transition. 

Another way to formulating this effective 5 dimen- 
sional theory is that, close to the minimum 8 = 0, if the 
scaling dimension of k z is fixed to be 1, then k x , k y ~ k8 
effectively have scaling dimension 2. Therefore expanded 
at the minimum the quadratic part of the free energy of 
$ reads: 



f , k l + 

F = J dk x dk y dk z ( — p 



7.2 



■r)\$ h 



(22) 



The total dimension is still 5, considering A^] = 
A [fey] = 2A[fe 2 ] = 2. All the other momentum depen- 
dent terms in the free energy are irrelevant. 

The symmetry of the lattice allows multiple degener- 
ate nodal points of function /(#, 4>) on the unit sphere 
labelled by solid angles. If the only nodal points are 
north and south poles 6 = 0, tt, which is allowed by 
the tetragonal symmetry of the lattice, the theory be- 
comes a precise 5 dimensional theory. However, the sym- 
metry of the system also allows four stable degenerate 
nodal points on the equator, for instance at (7r/2, mr/2) 
with n — ~ 3. Close to nodal points n = 0, 2, 
A[k z ] = A [fey] = 2A[fc x ] = 2, while close to nodal points 
n = 1, 3, A[k z ] = A[k x ] = 2A[k y ] = 2. Therefore 
the scattering between these nodal points complicates 
the naive counting of the scaling dimensions, although 
A[fc z ] = 2 is still valid. The transition in this case may 
still be a stable mean field transition, but more careful 
analysis of the loop diagrams is demanded to be certain. 
Let us denote the $ mode at (n/2, 0) and (ir/2, n) as $1 
and $*, and denote (ir/2,ir/2) and (tt/2, — 7r/2) modes 



F I ,//•■ c dk y dk z (^^ + k 2 x + r)|$ lifc | 2 



+ (^fi + k 2 + r)\^ k \ 2 



,k2 ^ a,k3 < * ) a,/€4 

a—1 i—1 
4 

i=l 
4 

i=l 
4 

S(Y J k l )g 2 ^ 1M ^ 1 . k2 ^> 1 . k3 ^2,k 4 - ( 23 ) 



The gi and g 2 terms describe the scattering between dif- 
ferent nodal points. To see whether the mean field tran- 
sition is stable, one can calculate the one-loop corrections 
to gi and g 2l and none of the loops introduces nonper- 
turbative divergence in the infrared limit. This analysis 
suggests that the Ising nematic transition is a mean field 
transition even with multiple nodal points of function 
f(8, </>) on the equator. 



V, SUMMARIES AND EXTENSIONS 

In this work we studied the global phase diagram of the 
magnetic order and lattice distortion of the Fe-pnictides 
superconductors. A two dimensional and three dimen- 
sional formalisms were used for 1111 and 122 materials 
respectively. The superconductivity was ignored so far in 
this material. If the quantum critical points discussed in 
this paper occur inside the superconducting phase, our 
results can be applied to the case when superconducting 
phase is suppressed. For instance, in 1111 materials, if a 
transverse magnetic field higher than H c2ja b is turned on, 
the field theory Eq. [2] and Eq. [6] become applicable. If 
the T c of the superconductor is lower than the ultraviolet 
cut-off of our field theory, the scaling behavior predicted 
in our work can be applied to the temperature between 
T c and the cut-off. Inside the superconducting phase, 
the nature of the transition may be changed. In 122 
materials, the ARPES measurements on single crystals 
indicate that the fermi pockets are fully gapped in the 
superconducting phase [29(, therefore the magnetic and 
nematic transitions are described by the z — 1, d = 3 
field theory Eq. [151 which is an extremely weak first 
order transition. In 1111 materials, although many ex- 
perimental facts support a fully gapped fermi surface, 
d— wave pairing with nodal points is still favored by the 
Andreev reflection measurements 27|, |28( . The nematic 
transition is the background of d— wave superconductor 
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is studied in Ref. [20J, [2J, |23j . 

In most recently discovered 11 materials 
Fei +y Se x Tei_ x , the SDW and lattice distortion are 



both different from the 1111 and 122 materials [24|. The 
SDW state breaks the reflection symmetries about both 
x = y line and x axis i.e. there are two different Ising 
symmetries broken in the SDW state, the ground state 
manifold is S 2 x Z2 x Z%. In this case, the classical 
and quantum phase diagrams are more interesting 
and richer, and since the order moments of the SDW 
in 11 materials are much larger than 1111 and 122 
materials (about 2fi B ), a lattice Heisenberg model with 
nearest neighbor, 2nd nearest neighbor, and 3rd nearest 
neighbor interactions ( J\ — J2 — J3) may be adequate in 
describing 11 materials, as was studied in Ref. [39j . 

Besides the quantum phase transitions studied in our 
current work, a quantum critical point is conjectured be- 
tween the P— based and As— based materials 35J, the 
field theory of this quantum critical point is analogous to 
Eq. [4] The formalism used in our work is also applicable 
to phase transitions in other strongly correlated mate- 
rials, for instance the spin-dimer material BaCuSi20g, 
which under strong magnetic field develops long range 
XY order interpreted as condensation of spin triplet com- 
ponent S z = — 1 [361 ]. This quantum critical point also 
has dynamical exponent z = 2, although the frequency 
linear term is from the Larmor precession induced by the 
magnetic field, instead of damping with particle-hole ex- 
citations. The frustration between the nearest neighbor 
layers in this material introduces an extra Ising symme- 
try between the even and odd layers besides the XY spin 
symmetry, therefore the quartic terms of this field theory 
are identical with Eq. [6] The RG equations of these quar- 
tic terms are much simpler than Eq. [5J because only the 
"ladder" like Feynman diagrams need to be taken into 
account [13] • We will study the material BaCuSi206 in 
detail in a future work [38]. 

We thank Subir Sachdev and Qimiao Si for heplful dis- 
cussions. We especially appreciate Bert Halperin for ed- 
ucating us about his early work on phase transitions on 
soft cubic lattice 13311. 
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